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^ ! Abstract 

We find a canonical form for pure states of a general multipartite 
Q_i| system, in which the constraints on the coordinates (with respect to 

a factorisable orthonormal basis) are simply that certain ones vanish 
and certain others are real. For identical particles they are invariant 
^ ' under permutations of the particles. As an application, we find the 

dimension of the generic local equivalence class. 



1 Introduction 

Recently considerable attention §, [| g @, |, [7|, |, |, [n|, pj, 0, |] has been 
devoted to the problem of describing the equivalence classes of states of a 
composite quantum system, where two states are regarded as equivalent if 
they are related by a unitary transformation which factorises into separate 
transformations on the component parts (a local unitary transformation). 
One approach to this problem is to specify a canonical form for states under 
local unitary transformations. For pure states of two-part systems, such a 
canonical form is given by the Schmidt decomposition 

\*) = Y t a i \<j> i )\il> i ) (1.1) 

i 

where the are a set of orthogonal states of the first subsystem, the \ipi) are 
a set of orthogonal states of the second subsystem, and the are positive real 
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numbers. The state |^) is thus expanded in terms of a factorisable basis of 
two-part states so that the number of non-zero coefficients is minimal. Acm 
et al [Q] have shown that there is an expansion of states of three qubits which 
has a similar property. In this note we will demonstrate such a decomposition 
for pure states of an n-part system, where the dimensions of the individual 
state spaces are finite but otherwise arbitrary. Then we use this to find the 
dimension of the generic local equivalence class. We will also comment on 
the relation of this decomposition to other proposed canonical forms. 



2 The generalised Schmidt decomposition 

We first state and prove the generalisation of the Schmidt decomposition for 
a multipartite state in which all the individual state spaces have the same 
dimension, since this is considerably simpler than the general case: 

Theorem 1. Let be a state vector in an n-fold tensor product space 
S\ <g) • • • ® S n where diniiSi = • • • = dimiS n = d > 2 and n > 3. Then for 
r = 1, . . . , n there is a basis ) : i = 1, . . . , d} of S r such that in the 

expansion 

l*>= E^l^)---!^) (2-1) 

the coefficients c^...^ have the following properties: 

2. Ci v ..i n is real and non-negative if at most one of the i r differs from d; 

3- |Cji...i| > l c ji...j,J — 3ri r = lj"-j71< 

Proof. Consider the real-valued function \ (^/\ (\4>^) ■ ■ ■ l^^))] 2 defined for 
unit vectors lying in the unit sphere S 2 ^ 1 in 7i r . As . . . , 

varies over the compact space S 2d ~ l x ■ ■ ■ x S 2d ~ l , this function attains a 

maximum at some point (\ip{^), ■ ■ ■ , • Let jlV^) : i = 1, . . . , dj be 

any orthonormal basis of Tt r containing and expand as in the state- 



ment of the theorem. Since 



(*i ( • • • • • • \4 n) ) 



2 

is stationary at 



\(f)^} = {ip^} for variations of \4>^) on the unit sphere, 



= <*| (|#> • • • l^-^l^)!^) • • • IVft) = for j > 1. 

(2.2) 
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Next, find the maximum of | (tf | (\<p^) • • • \<p {n) ))\ 2 as . . . , vary 

over unit vectors orthogonal to , . . . ,\ip[ n ^) respectively. Suppose the 
maximum occurs at ■ ■ ■ , IV^)) ■ Then, as before, in any expan- 

sion of |\&) in terms of orthonormal bases of the H r containing and 
\ip2^)i ( r = 1) • • • > n )' the coefficients will satisfy 

c 2 ... 2j 2...2 = for j > 2. (2.3) 

Continuing in this way, we define the basis vectors {tp^), • • • , IV^) f° r * = 
1, . . . , d— 1. Then the last basis vector \ip d ^) of Hi is determined up to phase. 

The reality conditions can be imposed as follows. Choose the phase of 
the basis element {tp^), i = 1, . . . , d — 1 so that 

wg(c d ...did...d) = (2.4) 

where the index % occurs in the rth place, and then fix Cd...d by choosing the 
phase of \ip d ^)- □ 

The general form of the theorem is rather more complicated than the 
above; to state it, we need to define the following sets of n-tuples. 

Let (Ji, . . . , Ijy) be the set of (n — l)-tuples . . . , i n _i) with 1 < i r < d r , 
excluding those of the form (i, . . . ,i) with 1 < i < d± and those of the form 
(di, . . . , d r , i, . . . , i) with d r < i < d r+ i and 1 < r < n — 2. We order 
these (n — l)-tuples in lexicographical order, so that In = (d±, . . . , d n -\). 
Let D = di ■ ■ ■ d n _i, so that N = D — d n _i + 1. We define A to be the 
set of n-tuples (ii, . . . , i n ) with . . . , i n -\) = Ik and i n = d n _\ + I where 
1 < k < min(iV, d n — d n -i) and k < I < d n — d n -\. 

We also define the following sets of n-tuples: 

B x ={(d 1 , d r -!,j, d r+1 , dn-i, d n -!) : 1 < j < d r , 1 < r < n - 2}, 
B 2 ={(di,...,d„_2,j, dn-i) : 1 < J < d„_ 2 }, 

#3 ={(^1, • • • , d n - 2 ,j, 1) : rfn-2 < j < rfn-l}, 

5 4 ={(2,2, . . . ,z) : 2<i<d 1 } 

U {(di, . . . , d r , i, . . . , i) : 1 < r < n — 1, d r < i < d r+ i} 
U {(di, . . . , d n _i, i) : d n _i < z < min(D, d n )}. 

Then A, Bi, . . . , B± are all disjoint. In terms of these sets, the general 
Schmidt decomposition can be stated as follows: 
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Theorem 2. Let |\&) 6e a state vector in an n-fold tensor product space 
Si <S> • • ■ <S> <S n where dimS r . = d r , with 2 < d± < ■ ■ ■ < d n , and n > 3. Then 
for r — 1, . . . , n there is a basis {IV^) '■ i = 1, ■ ■ ■ , d r } of S r such that in the 
expansion 

l*>= E^-l^)---!^) (2.5) 

the coefficients c^...^ have the following properties: 

1. Cn...ijii...i = if 1 < i < di and i < j. 

2. r ,/._.. = if d r < i < d r+ i and i < j, 1 < r < n - 2. 

3. cj = for every n-tuple index I in the set A. 

4. The coefficients with indices in the sets Bi {i — 1, . . . , 4) are real and 
non-negative. 

5. For i — 1, . . . , d n _i, define 

where r is such that d r < i < d r+ \. Then 

Proof. Consider the real-valued function |(^| • • • \4>^)) f defined for 

unit vectors \(f>^) lying in the unit sphere S 2dr ~ l in S r . As (|0^), ■ ■ ■ , \(f>^)) 
varies over the compact space 5' 2dl ~ 1 x • • • x S 2dn ~ l , this function attains a 
maximum at some point (l^i^)) • • • , IV'i™^))- Let {(V^) : i — 1, - ■ ■, d r } be 
any orthonormal basis of S r containing \tp^), an d expand as in the state- 



ment of the theorem. Since 



2 

is stationary at 



(^(l^V-l^) •••!#>) 

l ( / )<r ' > ) — IV^) f° r variations of \<fi^) on the unit sphere, 

= <*| • • • iVT^lVf >l^ +1) > • • • l^ n) >) = for j > I. 

(2.6) 

Next, find the maximum of | (^| • • • |^ (n) )) f as . . . , 

vary over unit vectors orthogonal to l^^}, ■ ■ ■ , l^i™' 1 ) respectively. Suppose 
the maximum occurs at ^|V4^)j • • • > IV^))- Then, as before, in any expan- 
sion of |\1>) in terms of orthonormal bases of the S r containing {ip^} and 
1^2^) ( r = 1) • • • > n )> the coefficients will satisfy 

c 2 ... 2j2 ... 2 = for j>2. (2.7) 
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Continuing in this way, we define the basis vectors IV^); ■ ■ ■ ■> IV^) f° r * = 
1, . . . , di — 1. Then the last basis vector (V'df ) of <Si is determined up to phase. 

Next, if the dimensions di are not all equal, we maximise (^/\ (jii^)^ 2 ' 1 ) • • • \4> <yn) ) 
as \4>^), • • • , \ (p^) vary orthogonally to the basis vectors already determined, 
to find basis vectors \tpf^), ■ ■ ■ , IV^); then we find ■ ■ ■ , \ip^) for % = 

di + 1, . . .,d 2 -l, and hence then |^ 3) ), . . . , |^- n) ) for i = d 2 , . . . , 4-1; 

and so on until we have \^ n ~^), for z = rf„_ 2 , . . . , — 1. The max- 

imisation at each step implies that the coefficients satisfy (2). 

Now, to fix the last d n — d n -\ + 1 basis elements of S n , we choose a set 
h — (ii, ■ ■ ■ , in-i) of n — 1 indices which is not of the form with 1 < 

i < d\ or (di, d 2 , ■ ■ ■ , d r , i, i, . . . , i) with d r < i < d r+ \ (so that Cj 1 ...i n _ 1 j has not 

yet been set to zero for any j), and maximise (\&| ^IV^) ' ' ' IV'ivTi^)!^ 7 ^)) 
with respect to vectors \<f>^) orthogonal to lip^), • • • , IV^-i-i)' thus finding 
iV'd^i)' then, choosing a different index set I 2 = (ji, ■ ■ ■ , jn-i), maximise 
(^1 [l^fi) ■ ■ - ^In-^)^)) 2 with res P ect to vectors |</> (n) ) orthogonal to 

■ ■ ■ , iV'd^-i)' anc ^ so on unt ^ we nave either exhausted the possible 
index sets h, I 2 , ■ ■ ■ , In or run out of space in which to vary the vector 
The coefficients will then satisfy (3). 

The reality conditions (4) can be imposed as follows. (By using a basis 
vector to fix a coefficient c we mean changing the phase of \tp) to make c 
real and non- negative.) First use \ip^) to fix c dl ... dn _ 1 x ; then use to 
fix c dld2 ... dn _ ldn _ 1 ; then use |^ r) ) (r = 1, . . . , n - 2; j = 1, . . . , d r - 1) to fix 
the coefficients in the set then use |'0 ? - ri ~ 1 ' > ) (j = 1, . . . , d n -i — 1) to fix 

the coefficients in S 2 and B 3 ; and finally use (j = 2, . . . , min(d n ,D) 

excluding j = c?„_i) to fix the remaining coefficients in i? 4 . □ 

This result can be expressed in terms of active transformations, with 
respect to fixed orthonormal bases {|#^)} of the state spaces S r , as follows. 
Two states \^ 2 ) in <Si <S> • • • <S> S n are said to be locally equivalent if 

= (tfi® ■■■®U n )\* 2 ) 
where U r is a unitary transformation acting on S r . Then we have 

Theorem 3. For n > 3, any state G S\ <8> • • • <g> <S n is locally equivalent to 

a state c h...i n \^) ' ' ' l^v^) where the coefficients c^...^ /iaue £/ie properties 
1-5 stated in Theorem 2. 
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As an example, we note the canonical form for a state of three qubits, 

a|000) + 6|011) + c|101) + d|110) + e|lll) (2.8) 

with b, c, d, e real. Acfn et al [[TJ have investigated tripartite states using 
different canonical forms which, like the canonical form proposed here, have 
five non-zero coefficients of which four are real. 

We note that for a state expressed in terms of a fixed orthonormal basis 

m (r) » as 

i^) = E^-iC ) )---i^: ) ), 

the coefficients of the basis elements = Yl^Jl^i^) defined in the first 

step of the above proof are the solutions of the generalised (nonlinear) "sin- 
gular value" equations 

^1 "'%t— lj^r+l ' 

where the Lagrange multiplier A is such that |A| 2 is the maximal value of 

K*| ■ •• I 2 - 

Let us now examine the dimension of the set of canonical forms and 
deduce the dimension of the generic local equivalence class. First consider 
the case where all the individual state spaces have equal dimension d. The 
number of zero coefficients in the canonical form, determined by condition 1 
of Theorem 1, is ^nd(d — 1) (one for each pair with i < j and for each 
position of j). The number of phases removed by condition 2 is n(d — 1) + 1. 
Hence the number of real parameters in the canonical form of Theorem 1 is 

2d n - nd{d - 1) - n(d - 1) - 1 = 2d n - n{d 2 - 1) - 1. 

We have not proved that states with different canonical forms are not lo- 
cally equivalent; it is conceivable that the number of parameters could be re- 
duced still further by local transformations. However, the difference between 
the number of parameters in the final canonical form and the dimension of 
the pure state space must be the dimension of the generic equivalence class, 
which is therefore at least n(d 2 — 1) + 1. But this is the dimension of the group 
of local unitary transformations, which can be identified with SU(d) n x U(l) 
if we collect together multiples of the identity in the final U(l). Since the 
local equivalence classes are orbits of this group, their dimension cannot be 
greater than the dimension of the group. Thus we have 

Corollary 1. If n > 3, the generic local equivalence class for a system of n 
d-state particles has dimension n(d 2 — 1) + 1. 
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Since a generic orbit has the same dimension as the group, the stabiliser 
of any state in the orbit has dimension zero. This gives 

Corollary 2. The generic pure state of a system of more than two equal-spin 
particles has a discrete stabiliser under the action of local unitary transfor- 
mations. 

The stabilisers of states of three qubits (d% = d 2 = d 3 = 2) were studied 
in |1J]. The second corollary generalises Theorem 1 of that paper. This treat- 
ment, however, gives no indication of which exceptional states have enlarged 
stabilisers. This question is being investigated in an alternative approach by 
one of us (HAC). 

For the general situation we must distinguish between the cases d n < D 
and d n > D where D = d\ . . . d n _i. In both cases, the number of zero 
coefficients imposed by conditions 1-2 of Theorem 2 is 

dl — 1 n n—2 d r +i — l « 

EEw-o+E E E«.-o 

8=1 s=l r=l i=d r s=r+l 

n n—2 n n—2 

=(di -l)) j d s + )Xd r+1 - d r ) 2j d s - nSt - }X n 

s=l r=l s=r+l r=l 

where S r 

n n-1 n— 1 

— — d s + d 2 r — S r — 5 n _x + d n d n ^i 

s=l r=l r=l 

n 

= \ d r (d r — 1) — ^5(6 + 1) where 5 = d n — d n -\. 

r=l 

If d n < D, the number of zero coefficients imposed by condition 3 is + 
while the number of phases removed by condition 4 is 

J2(dr - 1) + I- 

r=l 

Hence the number of real parameters in the above canonical form is 

n / n 

2]Jd r - J>^-1) + 1 

r=l \r=l 

which is the difference between the dimension of the state space and the 
dimension of the group G of local transformations. Thus in this case the 



- r)(S r+ i - S r ) 
= \d r (d r - 1) 
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dimension of the generic local equivalence class is the same as the dimension 
of G, as in Corollary 1. 

If d n > D, the number of zero coefficients imposed by condition 3 is 

TV 

J2(S — k + i) — N(S + 1) - \N{N + 1) 

k=l 

giving a total number of zero coefficients determined by conditions 1-3 as 

n 

\ d r (d r - 1) - \8{5 + 1) + \N{25 -N + l) 

r=l 

n 

=|^rf r (rf r -l)-lA(A-l) 

r=l 

where A = 5 — N +1 = d n — D. In this case there are no non-zero coefficients 
c ix ... in with i n > D, so the number of phases removed by condition 4 is reduced 
by A. Hence the total number of parameters removed, i.e. the dimension of 
the orbit, is at least 

n 

- 1) + 1 - A 2 = dim G - A 2 

r=l 

where G = SU^) x • ■ • x SU(d n ) x U(l). 

The fact that there are no non-zero coefficients c^...^ with i n > D means 
that the state is unaffected by unitary transformations of the nth particle 
which fix the first D basis vectors. Thus the stability group of the state 
contains at least this U (A) subgroup, and the dimension of the orbit cannot 
be greater than dimG — A 2 . It follows that the dimension of the orbit is 
exactly this, and the Lie algebra of the stability group is exactly that of 
U(A). 

Thus the general versions of Corollaries 1 and 2 are 

Corollary 3. In the general n-party system of Theorem 2, where n > 3, the 
generic orbit has dimension 

n 

^(d 2 -l) + l if d n <D, 

r=l 

n 

J2(dr-l) + l-(d n -D) 2 if d n >D, 

r=l 

where D = d\ - • • d n -i- 
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Corollary 4. In the n-party system of Theorem 2, a generic point has dis- 
crete stabiliser in the group of local unitary transformations if d n < D; oth- 
erwise the stabiliser is locally isomorphic to the unitary group U (d n — D) . 

3 Alternative generalisations of the Schmidt 
decomposition 

The canonical form of the previous section is noteworthy for the simplicity 
of the conditions on the coefficients in the expansion of the state vector in 
terms of a factorisable orthonormal basis: certain coefficients are zero, certain 
others are real. This is likely to make it most useful in practice. However, 
it is perhaps less theoretically appealing than some other generalisations of 
the Schmidt decomposition that have been proposed recently. For the sake 
of completeness, we review these alternatives here. 

For each constituent of the multipartite system, a basis of its individual 
state space is determined by its marginal density matrix: this is the basis 
defined by the conventional Schmidt decomposition of the state vector when 
the multipartite state space is regarded as a bipartite tensor product, one 
factor being the state space of the constituent being considered, the other 



being the tensor product of all the other state spaces. One of us []IT] and 
Brun and Cohen |2j have proposed that the tensor products of these one- 
particle states provide a natural basis for multipartite states. The resulting 
coefficients c^...^ satisfy 



11. ..I, — i,J r +i... l n 



for each r. 

Spekkens and Sipe |10| have suggested that a canonical state in each 



equivalence class could be taken to be that which minimises the Ingarden- 
Urbanik entropy 

Siv = ~Yl l c n-d 2l °gl c u-d 2 - 

ii—in 

They justify this as a generalisation of the Schmidt decomposition by showing 
that the IU entropy is minimised by the Schmidt normal form for bipartite 
states. For more than two constituent parts, however, little is known about 
these minima. 

To show that all three of these canonical forms are distinct, consider the 
tripartite state 



^) = ^=(3|000) + |011) + v / 2|lll)) 
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where each constituent system is a qubit (a two-state system) and, as usual, 
we label the basis states by the digits and 1 and abbreviate the product basis 
states as \abc) = \tpi }\ipl }\ipc )• This state is presented in the canonical 
form of section 2; not only does it satisfy the conditions of Theorem 1, but 
it can be shown (see Appendix) that it is obtained by the procedure of 
that theorem, i.e. the coefficient of 1 000) is maximal among states locally 
equivalent to (We note that there is a locally equivalent state with the 
coefficient of 1 000) equal to |; this satisfies conditions 1 and 2 of Theorem 1, 
but not condition 3.) By means of a transformation of the first qubit with 

matrix (^f^ vf+i))' ^ e s ^ a ^ e 1^) I s locally equivalent to 



m = 27! ( (v • " (V2 " 1)|100) 4 |m) 



which is in the canonical form of |TJJ and [0. Neither of these states minimises 
the IU entropy within their local equivalence class: under the infinitesimal 
local transformation in which the basis states of the first qubit transform by 
|0) i — ► |0) +e|l), |1) I— > |1) — e|0) and the states of the second and third qubits 
are kept fixed, the entropies of |^) and |$) change by 

6S m (V) = -jk el °Z 2 > 6S ™(®) = -elog(v^ + l). 

By minimising the IU entropy numerically we find in general that a few of 
the coefficients become much smaller than the others without vanishing 
exactly. 
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CLVqWq + bV\W\ = 


= Am 


CV\Wi = 


= Xui 


HUqWo = 


= Xvo 


bu Wi + cuiWi = 


= Xvi 


au v = 


= Xw^ 


TwqVx + CU\V\ = 


= XvJi 



Appendix 

We will show that for a state 

= o|000) + 6|011) + c|lll) 

with \a\ > 4j, a is the maximal value of the coefficient of 1 000) among states 
locally equivalent to |\&). 

Proof. The equations (|Z9|) for a stationary value of | ($\ |</> (2) ) |</> (3) )) | 2 

become 

(3.1) 
(3.2) 
(3.3) 
(3.4) 
(3.5) 
(3.6) 

where l^ 1 )) = u \0) + |0 (2) ) = v o \0)+v 1 \l) and |0( 3 )) = w \0) + w^l). 

Clearly there is a solution A = a, u = v = w = (1,0). We have to show that 
any other solution has |A| 2 < \a\ 2 . Using ( |3.5| ) and ( |3.6| ) to eliminate w and 
Wi, eqs. (|3.1|) and (|3.2|) become 

(|a| 2 K| 2 + |6| 2 M 2 -|A| 2 ) «o + ScNV = 0, 
bc\v 1 \ 2 u,+ {\c\ 2 \v 1 \ 2 -\X\ 2 ) u 1 = 0. 

For (u Q ,Ui) 7^ (0,0), it follows that 

F(|A| 2 ) = |A| 4 - |A| 2 (H 2 + (1 - 2|a| 2 )h| 2 ) + |a| 2 |c| 2 h| 2 M 2 = 

since |a| 2 + \b\ 2 + \c\ 2 = \v \ 2 + |wi| 2 = 1. Now 

F(|a| 2 ) = |a| 2 |t; 1 | 2 (2|a| 2 -l + |c| 2 |t;o| 2 ). 

which is positive if \a\ > -j= (unless V\ = 0), and the gradient of the quadratic 
F at |A| 2 = \a\ 2 is also positive. It follows that the zeros of F, and therefore 
any stationary values of | {\4>^)\4>^)\4>^)) | other than |a| 2 , are less than 
\a\ 2 . □ 
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